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Abstract 

The n-th moments of the nuclear structure function F^(x,Q 2 ) are an- 
alyzed using the off-shell kinematics appropriate to describe deep inelastic 
electron-nucleus scattering within the impulse approximation. It is shown 
that off-shell effects are sensitive to the form of both the nuclear spectral 
function and the nucleon structure function F^ \x,Q 2 ), and can be as large 

as ~ 10% at Q 2 ~ 2 (GeV/c) 2 . 
PACS numbers: 13.60Le, 25.30Fj, 25.30Rw 
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Nuclear effects in deep inelastic lepton-nucleus scattering (DIS) have been extensively 
studied, both experimentally and theoretically, over the last two decades (for a review see 
refs. [0,0)- The body of available data provides clearcut evidence that the nucleus cannot be 
simply described as a collection of nucleons on the mass shell. Proton knock-out experiments 
in which the outgoing particle is detected in coincidence with the scattered electron have 
clearly shown that the nuclear spectral function P(p, E), yielding the probability of finding a 
nucleon with momentum p and removal energy E in the target nucleus, extends to very large 
values of |p| (> 800 MeV/c) and E (> 200 MeV) ||. Hence, the theoretical description 
of lepton scattering off fast off-shell nucleons must be regarded as a prerequisite for the 
understanding of nuclear DIS data. 

A pioneering study of the electromagnetic response of a bound nucleon, focused on the 
intermediate energy domain, has been carried out in the 80's by de Forest ||]. More recent 
theoretical analyses of off-shell effects in nuclear DIS 0|| have been stimulated by the 
availability of new experimental information on the Q 2_ dependence of the moments of the 
nuclear structure function F^(x,Q 2 ), showing sizable deviations from the prediction of the 
simple convolution model ||. 

In this paper we study the relevance of nuclear effects to the Q 2 -dependence of the nuclear 
structure function F^(x,Q 2 ) and its moments using the appropriate off-shell kinematics 
and a realistic nuclear spectral function, resulting from ab initio microscopic many-body 
calculations 0. 

Let us consider electron scattering off a nucleon carrying four momentum p = (p , p) with 
p 2 ^ m 2 , m being the nucleon mass. As it is well known (see, e.g., ref. ||), the differential 
cross section for the inclusive process e + N — > e' + X, in which the hadronic final state 
is undetected, can be written in the form 



dVldE e , F(k,p) Q 4 E e 

where a 2 = e 2 /4ir = 1/137, q 2 = v 2 — |q| 2 = — Q 2 is the squared four-momentum tran- 
fer and E e and E e > are the initial and final electron energy, respectively. The function 



F(k,p) = v re i = \(k/E e ) — (p/po)|, k = (.E^k) being the four momentum carried by 
the incoming electron, describes the incident flux. In the case of ultrarelativistic electrons 
it takes the simple form F(k,p) = (kp)/(poE e ). Note that for an isolated stationary nucleon 
F = 1. 

The electromagnetic structure of the electron and the nucleon is described by the ten- 
sors L^ u and W^, respectively. Expanding the hadronic tensor in terms of the available 
relativistic invariants and taking into account current conservation, the contraction L^W^ 
can be readily written in terms of two nucleon structure functions W\ and W 2 : 



ITW^ = 2\2 W x {kk') + \2(kp)(k'p) - p 2 (kk' 



(2) 



where k! = k + q denotes the four momentum carried by the outgoing electron and (kk') = 
2E e E e i sin 2 (0/2), 6 being the electron scattering angle. In the case of scattering off a free 
nucleon at rest p 2 = m 2 , {kp)(k'p) = m 2 E e E e r and substitution of eq.@ into eq.([l|) yields 
that standard result 



d 2 a N 
dVtdE P i 



cm 



a 

W 2 + 2W X tan 2 - 



(3) 



where <jm = « 2 cos 2 (9 / 2) / [4E 2 sin 4 (6'/2)] is the Mott cross section and the nucleon structure 
functions depend upon q 2 and (qp) only. In general, the inclusive electron-nucleon cross- 
section contains additional terms and can be cast in the form 



d 2 (iN 
dQdE P 



F(k,p) 



W 2 {1 - A) + 2 W^tan 2 - 



(4) 



with 



A(k,k',p) 



(kp) (k'p) 



mEf, mEfi 



cos 



P 
m 



- tan - 
2 / 2 



(5) 



and the structure functions W\ and W 2 depend upon the three relativistic invariants q 2 , (qp) 



and p 2 . For an isolated stationary nucleon p 2 
is recovered. 



m 



F(k,p) = 1, A(k,k',p) = and eq.(|) 
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The electron-nucleon cross section of eq. (|4]) can be used to obtain the inclusive electron 
nucleus cross section within the impulse approximation from 0: 

da a 



/ d 2 a N 
\dQdE e . 



(6) 



dVtdE e i J \po 

where P(p) = P(\p\,m — p Q ) denotes the nuclear spectral function. The interpretation 
of the above equation is straightforward. The electron-nucleus cross section is written as 
the incoherent sum of the cross sections associated with scattering off individual nucleons, 
weighted with the nuclear spectral function and stripped of the flux factors, since the incident 
flux has to be defined with respect to the target, i.e. to the nucleus as a whole. Note that 
the nuclear cross section of eq.@ can be recast in a form similar to eq.(|3]), with the nucleon 
structure functions W\ and W 2 replaced by two nuclear structure functions W A and W A . 

Due to the fact that P(p) is a steeply decreasing function of |p|, the main contribution 
to the |p| integration involved in eq.(^) comes from the region of p 2 small compared to m 2 , 
where the use of a nonrelativistic nuclear spectral function is quite reasonable. Furthermore, 
as a first approximation, the functions W\ and W 2 , whose p 2 -dependence is unknown, can 
be estimated assuming Wi t2 (q 2 , {qp),p 2 ) — Wi t2 (q 2 , {qp),™ 2 ), and Wi^iq 2 , (qp),m 2 ) can be 
extracted from electron-proton and electron-deuteron scattering data. 

Using eqs.([|) and (j^j the nuclear structure function F 2 A = vW 2 A can be related to the 
nucleon structure function F 2 = vW 2 through 

[■A 



F 2 A (x, Q 2 ) = J dz <h(z, Q 2 )F 2 g, Q 2 ) , (7) 



where 



Uz,Q 2 )=zJd^p {^)P{P) [1-A(M',P)] 5 { z -^^) . (8) 

Ma being the mass of the target nucleus. 

Note that setting A = eq.(^) reduces to the standard convolution formula (see, e.g., 
ref [0): 

F 2 A (x,Q 2 ) = jy z f A (z,Q 2 )F 2 g,Q 2 ) , (9) 



with 

Let us now focus on the calculation of the n-th moment of the nuclear structure function, 
defined as 

M^(Q 2 ) = [ A dx x n ~ 2 F 2 A (x,Q 2 ) . (11) 

Using the convolution model result of eq.@ and assuming the validity of the Bjorken limit, 
which amounts to disregard the Q 2 -dependence of f A (z, Q 2 ) of eq. (pT0|) , eq. (|TlD can be rewrit- 
ten in the simplified form 

M$ (Q 2 ) = j* dz £ dy z n - 1 f A (z)y n - 2 F 2 (y, Q 2 ) , (12) 

leading to the factorized expression 

M^Q 2 ) = f^M^Q 2 ) . (13) 

In the above equation M^(Q 2 ) is the n-th nucleon moment and 

f { A +1) = I* dz z"- 1 f A (z) , (14) 

J 



whereas f A (z) denotes the Q 2 — > oo limit of Ja(z,Q 2 ) [ T0|1 . 

Substitution of the nuclear structure function given by eqs. (|7]) and @ into eq. (|TT|) shows 
that the factorization property exhibited by the n-th moment evaluated within the convo- 
lution model breaks down when the Fermi motion and binding of the struck nucleon are 
properly taken into acount. 

The breaking of factorization can be best observed studying the ratio between the nuclear 
and nucleon n-th moments, defined as 

According to the convolution model, R n (Q 2 ) is Q 2 -independent, hence the study of its Q 2 
dependence may provide useful information on the relevance of the factorization breaking 



terms arising from the inclusion of off-shell kinematics for the struck nucleon in the calcula- 
tion of F 2 A from eqs.(0)-(|). 

The results of our numerical investigation of the effects of factorization breaking on the 
ratio R 5 (Q 2 ) are shown in fig.|I]. The solid curve has been obtained using F 2 (x, Q 2 ) given by 
eqs. (0) and (§) and the nuclear matter spectral function of ref . |7[] , whereas the data points, 
taken from ref. ||, have been extracted from the analysis of several experiments carried 
out at CERN ]ll| and SLAC |]l2| , |l3f using 56 Fe targets. It appears that the calculated 
Rb{Q 2 ) decreases in the range 2 < Q 2 < 10 (GeV/c) 2 , while becoming almost constant at 
Q 2 > 10 (GeV/c) 2 . The maximum deviation from the prediction of the convolution model, 
shown by the horizontal line in fig. [L], is ~ 15%. Our results appreciably depend upon the 
form of the nucleon structure function F 2 (x,Q 2 ), particularly at Q 2 < 10 (GeV/c) 2 . In this 
region we have employed the parametrization proposed in refs. |14],|15j , generally referred to 
as CKMT model, that provides a consistent formulation of the structure function at low 
Q 2 and gives a good description of the recent data from HERA |p!6|-|20|]. According to the 
CKMT model, based on Regge theory, the nucleon structure function F 2 can be split into 
two parts: the singlet term, corresponding to the Pomeron contribution, and the nonsinglet 
term, corresponding to the secondary Reggeon. Both terms have the nonpertubative Q 2 - 
dependence. At Q 2 > 10 (GeV/c) 2 we have used the so called GRV model, proposed in ref. 
nj, which reproduces the Q 2 -evolution of F 2 (x,Q 2 ) quite well. 
We have also used our approach to calculate the quantity 

,2M -(V<*«) 



T n(Q ) 



(16) 



where d n = ■y n /2j3 , j n is the so called non singlet (NS) anomalous dimension and 
(3q =11 — (2/3)nf, rif being the number of flavors. 

The solid line of fig. ^| shows t^(Q 2 ), obtained from eqs.(0)-(|) and (|TT|), for Ql = 12.5 
(GeV/c) 2 . It exhibits a deviation from the log(Q 2 /A 2 ) behavior predicted by leading order 
perturbative QCD and the factorization theorem. This deviation is a consequence of both 
higher twist corrrections, included in the CKMT parametrization of the nucleon structure 

6 



function, and nuclear effects. In order to illustrate the relative importance of the two sources 
of non logarithmic Q 2 -dependence of t^(Q 2 ), we also show, by the dashed line, the results 
obtained setting A = in eq.(§). Comparison between the solid and dashed lines shows 
that nuclear effects reach a maximum of ~ 10 % at Q 2 = 2(GeV/c) 2 . 

The results presented in figs. [I] and |2] can be related to each other noting that, defining 
the nucleon analog of t£(Q 2 ) as tJ?(Q 2 ) = [Mf (Q 2 ) / 'M$ ' (Q 2 )]- 1/d5 , we can write 



ri{Q 2 ) 



R,{Q 2 



l/ds 



(17) 



Numerical calculations show that t^(Q 2 ) is very close to the quantity represented by the 
dashed line of fig. |2|. Hence, to a very good approximation, the above equation provides 
a relationship between R 5 (Q 2 ) of fig. [I] and the ratio of the quantities represented by the 
dashed and solid lines of fig. 0. Since the experimental errors in t^(Q 2 ) are smaller than 
those associated with Rn(Q 2 ) fig. H provides a better ilustration of our results. Note that 
we obtain similar results for all values of n (n=3-8) for which data is available f5|. 

The Q 2 -dependence of n-th moment ratios and the relevance of off-shell corrections has 
been recently discussed by Cothran et al [fj]. In ref. || the structure function of an on-shell 
nucleon, F° N is written in terms of the relativistic invariant vertex function <&(S) as 



?on. 



x) 



X 



d 2 p t 



x 2 (S — m 2 ) 



(18) 



l-xJ 2(2tt) 3 

where S = p 2 /(x(l — x)) + A 2 /(l — x) + fi 2 /x, /i and A being the masses of the struck 
and spectator quarks respectively Applying Old- Fashioned- Perturbation-Theory to a 

system consisting of a constituent of mass fi and unit charge plus a neutral core of mass A 



one finds [22 



F 2 {x) 



x 



d 2 p t 



d\ 2 



xS 



-. 2 



(19) 



l-x J 2(2vr) 3 

where S = m 2 — S and Q\(xS)/S can be identified with the Fock space wave function in 
the infinite momentum frame. Comparison to eq.(|T9"D shows that eq.(^) can be obtained 
carrying out an integration over the squared mass of the spectator system, A 2 , and is strictly 



applicable only in the x — > 1 limit (see, e.g., ref. ||). Our numerical results show that the 
integrations involved in the calculations of the nuclear structure function (see eq.(|7])) receive 
non negligible contributions from the region of small (x/z) (for example, at Q 2 = 2 (GeV/c) 2 
and x=0.7 more than 50 % of the integrated strength comes from the region (x/z) < 0.7). 
Hence, the results of the approach of ref. |5[] should be taken with some caution. 

In conclusion, we have shown that the use of off-shell kinematics, appropriate to describe 
lepton scattering off bound nucleons, in the calculation of the nuclear structure function 
F 2 (x, Q 2 ) leads to a sizable breakdown of factorization. The amount of the effect depends 
upon both the form of nuclear spectral function and the model of the nucleon structure 
function F 2 (x, Q 2 ). 

Our approach is based on a standard nonrelativistic many-body treatment of nuclear 
dynamics, involving no adjustable parameters. Pionic and relativistic corrections may also 
play a role in this context, but unfortunately our formalism has not yet been developed to 
include their contributions in a consistent fashion. 

We find that the the maximum deviation associated with off-shell effects is ~ 10% and 
occurs at Q 2 ~ 2 (GeV/c) 2 , while at Q 2 > 5 (GeV/c) 2 the effect becomes negligibly small. 
Our calculations, based on a realistic spectral function and the CKMT model for the nucleon 
structure function, suggest that nuclear effects may be larger than previoulsy estimated || , 
and that the extraction of higher twist corrections to F 2 (x,Q 2 ) from nuclear data may be 
questionable. 

The authors gratefully acknowledge very helpful discussions with W. Weise. 



8 



REFERENCES 

t Supported by the Russian Foundation for Fundamental Research under grant N 96-02- 
16126 

[1] M. Arneodo, Phys.Rep. 240 (1994) 301. 

[2] D.F. Geesaman, K. Saito and A.W. Thomas, Ann. Rev. Nucl. Part. Sci. 45 (1195) 337. 

[3] O. Benhar, V.R. Pandharipande and S.C. Pieper, Rev. Mod. Phys. 65 (1993) 817. 

[4] T. de Forest, Jr., Nucl.Phys. A392 (1983) 232. 

[5] CD. Cothran, D.B. Day and S. Liuti, Phys.Lett. B421 (1998) 46. 

[6] A.S. Rinat and M.F. Taraign, preprint WIS-99/13 (1999) 

[7] O. Benhar, A. Fabrocini and S. Fantoni, Nucl. Phys. A505 (1989) 267. 

[8] F.E. Close, An Introduction to quarks and partons (Academic Press, New York, 1979). 

[9] O. Benhar, S. Fantoni, G.I. Lykasov and N.V. Slavin, Phys.Rev. C55 (1997) 244. 
[10] O. Benhar, S. Fantoni, G.I. Lykasov and N.V. Slavin, Phys.Rev. C57 (1998) 1532. 
[11] J.J. Aubert et al, Phys. Lett. 123B (1983) 275. 
[12] J. Gomez et al, Phys. Rev. D 49 (1994) 4348. 
[13] J. Arrington et al, Phys. Rev. C 53 (1996) 2248. 

[14] A. Capella, A.B. Kaidalov, C. Merino and J. Tran Thanh Van, Phys.Lett. B337 (1994) 
358. 

[15] A.B. Kaidalov and C. Merino, Eur. J. Phys. CIO (1999) 153. 
[16] NMC, P. Amaudruz et al, Phys. Lett. B259, (1992) 159; 
[17] HI Collaboration, I. Abt et al, Nucl. Phys. B407 (1993) 515. 

9 



[18] ZEUS Collaboration, M. Derick et al, Phys. Lett. B316 (1993) 412. 

[19] HI Collaboration, C. Adloff et al, Nucl. Phys. B497 (1997) 3. 

[20] ZEUS Collaboration, J. Breitweg et al, Phys. Lett. B407 (1997) 432. 

[21] M. Gliick, E. Reya and A. Vogt, Z. Phys. C67 (1995) 431. 

[22] D. Sivers, S.J. Brodsky and R. Blankenbecler, Phys.Rep. 23C (1976) 1. 



10 



FIGURES 

FIG. 1. Q 2 -dependence of R 5 (Q 2 ) = M 5 A (Q 2 )/Mf (Q 2 ). The solid curve shows the results 
of our approach, whereas the horizontal line represents the prediction of the convolution model 
(eqs.(]9|) and (|i~3|)). The data is taken from ref. ||. 

FIG. 2. Q 2 -dependence of [M^(Q 2 )/M^(Q 2 ) )}- 1 / d5 . The solid line corresponds to the full 
calculation, whereas the dashed line has been obtained neglecting off-shell corrections, i.e. setting 
A = in eq.(§). 
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